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The phase diagram of one-dimensional earth-alkaline fermionic atoms and ytterbium 
171 atoms is investigated by means of a low-energy approach and density-matrix 
renormalization group calculations. For incommensurate filling, four gapless phases with 
a spin gap are found and consist of two superconducting instabilities and two coexisting 
■ bond and charge density-waves instabilities. In the half-filled case, seven Mott-insulating 

phases arise with the emergence of four non-degenerate phases with exotic hidden 
orderings. 
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The recent experimental progresses achieved in trapped ultracold atomic gases 
augur great opportunities to explore the physics of strong correlations in clean 



5-H ' systems and within a wide range of parameters, owing to optical lattices and the 

high tunability of Feshbach resonances. The attention of the cold atoms community 
has recently been called to fermionic earth-alkaline atoms such as strontium (Sr) 
and to atoms endowed with a similar electron structure, such as y tterbium (Yb), 
which have been cooled down to reach the quantum degeneracy 'l^. The interest in 
those atoms stems from the presence of a ground state ^Sq ("5") and a metastable 
excited state ^Pq ("e") between which transitions are forbidden; moreover, both 
states have zero electronic angular momentum J, so that the nuclear spin / is 
decoupled from the electronic spin 131. The scattering lengths are independent of the 
nuclear spins and this results in systems with an extended SU(Af = 2/+1) symmetry 
that can be realized without any fine-tunning; N can be as large as 10 with ®''Sr 
(/ = 9/2). In this respect, earth-alkaline cold atoms seem to be very promising 
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for the experimental realiz atioii of high-symmetry systems, and the investigation of 
exotic many-body physics | 4 | 5 | 6 | 7 | 

In this paper, we will focus on the one-dimensional (ID) N = 2 case, that is 
on systems of earth-alkaline like atoms with a nuclear spin 1 — 1/2 such as -'^'^^Yb 
atoms or earth-alkaline atoms where only two nuclear spin states have been trapped. 
The general N case will be investigated elsewhere. Besides, in addition to the U(l)c 
charge symmetry due to the conservation of the total number of atoms, there is 
also a U(l)o orbital symmetry, as a result of the conservation of the number of 
atoms in each state {g and e). These systems display thus an SU(2)sXU(l)cxU(l)o 
continuous symmetry. Moreover, for the sake of simplicity, we assume an extra 
orbital Z2 symmetry between the two orbital levels e and g so that they play a 
similar role. In that case, the effective Hamiltonian model describing such systems 
loaded in ID optical lattice is: ^ 

i,la i,l i.l 

+ ^ H + ^cx 4a,^cI/^,^Cg/3,^Cea,^, (1) 

where a,/? =t7i is the SU(2) spin index and I = g,e represents the two orbital 
levels. In Eq. ([1]), ■ are the four-components fermionic creation operators on 

site i, and ni^i — J2a "^L i'^ia i denotes the density operator on orbital level 
I = e,g and site i. We intend to investigate the zero-temperature phase diagram of 
model ([T]) for incommensurate filling and at half-filling. For this purpose, we will 
follow a weak-coupling approach, supported by numerical simulations for moderate 
couplings, using the Density Matrix Renormalization Group (DMRG) algorithmic]. 



1. Incommensurate filling 

Since we are interested in the low-energy properties of model ([1]), the starting point 
of our approach lies in the weak-coupling regime \U, V, Vex\ <C t. In this regime, only 
the long-distance properties of the model matter and we perform a continuum limit 
by linearizing the dispersion relation for non-interacting fermions about the two 
Fermi points {±kp). At incommensurate fillings, there is a spin-charge separation: 
the charge and the spin, orbital degrees of freedom are decoupled and can be treated 
separately. The charge degrees of freedom are gapless and exhibit the characteristics 
of the Luttinger liquid The low-energy properties of the remaining spin and 
orbital degrees of freedom can be described in terms of six Majorana (real) fermions 
^RL — Ij • • • 76) as in the spin-orbital ladder The interacting part of the 
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corresponding effective Hamiltonian density reads as follows: 
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with gi = -ao(;7 + Vcx), 32 = -aoV', 53 = -ao(C/- V^x), 34 = -ao(T^ - 2V;x), and 
35 = ao(C^ — 21^ + Vox), ao being the lattice spacing. The one-loop renormalization 
group (RG) equations for that model can be found in Ref. |TT] The RG analysis 
leads to the emergence of four different spin-gapped phases. Two phases have 
coexisting bond and density wave instabilities and can be distinguished by their 
properties under the Z2 exchange symmetry (e O 17). Spin-Peierls (SP) and 
charge-density wave (CDW) instabilities are even under Z2 while SP^r and CDW^r - 
or orbital-density wave (ODW) - instabilities are odd. Their order parameters are: 
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l,a 



^ g-2.fc,ao t 



La 



SP _ g-2ifcFao 
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e-^^''-'^°Y.'lcl,^^la,+l. (3) 

l,a 



with ei — ±1 for / = g,e. The SP instability has dimerised pairs of atoms in the g 
or in the e states; the CDW instability has alternating fully occupied and empty 
states, and the CDW^r (or ODW) instability is an alternation of sites with two 
atoms in the g state and sites with two atoms in the e state. Figure [3] (next section) 
shows a pictural representation of these instabilities. 

The two other phases are superconducting instabilities and differ in the Z2 
exchange symmetry of their pairing operators: 



/-/iBCSs _ 
QBCSd 



Cgt,i+lCeUi - Cg4„i+lCet,i - (e O g), 
- Cgf^iCe_l.i — Cg^^iCg-f^i. (4) 

The zero-temperature phase diagram is sketched in Fig. [1] depending on the sign of 



2. Half-filling case 

2.1. Weak- coupling approach 

At half-filling, there is no spin-charge separation anymore since umklapp processes 
couple these degrees of freedom in sharp contrast to the half-filled Hubbard 
chain The interplay between charge and spin-orbital degrees of freedom can 
be investigated by introducing two additional Majorana fermions ^ for the 
charge sector. The umklapp processes generate new terms which add to the effective 
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Fig. 1. Phase diagram of model JTJ at incommensurate filling for Vex = — O.Olt (left) and Vej 
O.Oli (right). 
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where = ao{U +2V-Vc^), gr = aaiV-Vc^), gg = aoU, and gg = ao{V+Vc,c). The 
anal ysis of the one-loop RG equations for the Majorana models (OE]) is detailed in 
Ref. L!^. The resulting zero-temperature phase diagram is depicted in Fig. [51 
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Fig. 2. Phase diagram of model ([TJ at half-filling for Vex = -0.03i (left) and Vex = O.Oli (right). 
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Seven Mott-insulating phases have been found in the phase diagram, aU fully 
gapped in charge, orbital, and spin degrees of freedom; three of them are two-fold 
degenerate and four are non-degenerate. As regards the three degenerate phases, 
depicted in Fig. [31 their order parameters were given in Eq. one is the SP 
phase (Of^), a second one is the CDW phase {Of^"^), and the third one is the 
ODW phase (Ofow ^ ^cdw^^^ incommensurate case, these orders have a 

power-law decay while true long-range ordering occurs in the half-filled case. 




Fig. 3. (color online) The blue atoms are in the ground state, and the yellow color represents the 
excited state, (a) The SP phase contains dimerised pairs of atoms in the g or in the e states; (b) 
the CDW phase swings between fully occupied and empty states; (c) the ODW is an alternation 
of sites with two atoms in the g state and sites with two atoms in the e state. 



The first non-degenerate phase is the rung triplet (RT) phase, so called since it 
is related to the RT phase of the two-leg spin-1/2 ladder, with the orbital states 
g, e playing the role of the legs. This gapped phase is adiabatically connected to the 
Haldane phase of the spin-1 Heisenberg chain 1^ and displays a dilute order which 
consists in an alternation of S*^ ^ Sg + — +1 (with S'f — [cj^ ^c^^ i ~ '^li i^ii il/"^) 
and — —1 states with an arbitrary number of S*^ = states in between (see 
Fig. H (a)). The second non-degenerate phase is a rung singlet (RS) phase where 
the two orbital states form a spin singlet state on each site (Fig. |4] (b)). The third 
phase is called Haldane charge (HC) phase since it belongs to the same type as the 
Haldane (spin) phase above, but now, the correct quantum number to characterize 
the dilute order is the occupation number rii. States with rii = A and = 
alternate, and there is an arbitrary number of singlet (n^ = 2) states in between. It 
is easily seen (Fig.|3](c)) that HC interpolates between the CDW and the RS phase. 
This HC phase has been first identified in the context of ID spin-3/2 cold alkaline 
fermions^!^. Finally, the last non-degenerate phase is called Haldane orbital (HO) 
and the alternation takes place between states with (n^ — ne)/2 = ±1 with again 
an arbitrary number of spin singlets {ug — ne)/2 = in between (Fig. 111(d)). This 
phase interpolates between the ODW and the RS phases. 
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Fig. 4. (color online) The blue atoms are in the ground state, and the yellow color represents the 
excited state, (a): The RT phase displays alternating = ±1 states with an arbitrary number 
of = (spin triplet) states in between; (b) the RS phase consists in spin singlets formed on 
each site by two atoms, one in the g state and the other in the e state, (c) the HC phase shows 
alternating rii = 4, states with an arbitrary number of singlets in between; (d) the HO phase 
exhibits alternating (rig ~ne)/2 = ±1 states with again an arbitrary number of singlets in between. 



2.2. Numerical approach 

We now carry out numerical calculations using DMRG in order to investigate 
the various phase diagrams, and we focus on the half-filled case. For a finite 
system of size L, we can fix several quantum numbers: the total number of 
particles Qc — ; ~ ^L, the z-component of the total spin of the whole system 
'5'tot = J^ii'^ei + i): a-s well as the population difference between the two orbital 
levels Qo = '^i{n-e,i ~ n-g,i)- Typically, we keep up to 1600 states, which allow to 
have an error below 10~^, and we use open boundary conditions (OBC). 

The low-energy approach predicts seven insulating phases in the phase diagram. 
In order to identify the degenerate Mott phases (CDW, ODW, and SP), we 
compute local quantities (bond kinetic energy, local density, etc.). In order to detect 
non-degenerate Mott phases, we investigate the presence or absence of various 
edge states with quantum numbers {Qc, S^^^,Qo), respectively equal to (L, 1,0), 
(1 + 2, 0, 0), and (L, 0, 1) for the RT, HC, and HO phases. The RS phase is identified 
as being non-degenerate and does not exhibit any kind of edge states with OBC. 

The resulting zero-temperature phase diagrams are shown in Fig. [5] for two 
choices of Vex- The overall topology of these two phase diagrams is in agreement 
with the low-energy results plotted in Fig. [5] since we confirm the presence of the 
seven insulating phases. There is a fair concordance for positive V^x- However, there 
is a little discrepancy between the two approaches for negative Vex- the small SP, HC 
and HO pockets do not match exactly. Still, we emphasize that the weak coupling 
phase diagram was obtained for much smaller values of | Vex \ ', furthermore, the weak 
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coupling approach shows that the SP pocket shifts and then quickly disappears for 
increasing |Vox|, in favor of the HC and HO phases. This indicates that the position 
of the SP pocket is not fully significant and strongly depends on the value of Vox- 
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Fig. 5. Numerical phase diagram of model JTJ at half-filling for Vex = —t (left) and Vex = t (right). 
We fix t = 1 as unit of energy. 



3. Concluding remarks 

We have investigated the phase diagram of ultracold ID fermionic systems of 
earth- alkaline like atoms with nuclear spin (/ = 1/2), at incommensurate filling 
and at half-filling. We showed that at incommensurate filling, while the charge 
degrees of freedoms remain gapless, the spin and orbital degrees of freedom are all 
gapped. The phase diagram displays four different phases with coexisting SP and 
CDW or with BCS instabilities. 

At half-filling, seven Mott-insulating phases arise. Three of them are two-fold 
degenerate (SP, CDW, and ODW phases) and four are non-degenerate (RT, RS, 
HC, and HO phases). Among the non-degenerate phases, three are of the Haldane 
type (related to the physics of spin chains) and show hidden dilute order. We hope 
that, in the light of the recent experiments of earth-alkaline-atoms and ^^^Yb atoms, 
future experiments will disclose part of the richness of the phase diagram presented 
in this paper. 
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